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Abstract 

Compressible, subsonic flow through a diffusing S- 
duct has been experimentally investigated. Benchmark 
aerodynamic data are presented for flow through a repre- 
sentative S-duct configuration. The collected data would 
be beneficial to aircraft inlet designers and is suitable for 
the validation of computational codes. Measurements of 
the three dimensional velocity field and total and static 
pressures were obtained at five cross-sectional planes. 
Surface static pressures and flow visualization also helped 
to reveal flow field characteristics. All reported tests were 
conducted with an inlet centerline Mach number of 0.6 
and a Reynolds number, based on the inlet centerline ve- 
locity and duct inlet diameter, of 2.6 x 10 6 . The results 
show that a large region of streamwise flow separation oc- 
curred within the ducL Details about the separated flow 
region, including mechanisms which drive this compli- 
cated flow phenomenon, are discussed. Transverse veloc- 
ity components indicate that the duct curvature induces 
strong pressure driven secondary flows, which evolve into 
a large pair of counter-rotating vortices. These vortices 
convect the low momentum fluid of the boundary layer 
towards the center of the duct, degrading both the unifor- 
mity and magnitude of the total pressure profile. 
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normalized Mach number 
static pressure 
total pressure 
center line arc length 
duct cross-sectional radius 
duct curvature radius 
Reynolds number 
law-of-the-wall ordinate 
cartesian coordinate system 
law-of-the-wall abscissa 
boundary layer thickness 
displacement thickness 
momentum thickness 
duct substended angle 
cross-stream polar angle 


Subscripts 


centerline 

station at beginning of curvature 
station at end of curvature 
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duct cross-sectional area 

sonic velocity 

static pressure coefficient 

total pressure coefficient 

inlet duct diameter 

boundary layer shape factor ( 61 / 62 ) 

Mach number 
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Introduction 

A IRCRAFT propulsion systems often use diffusing S- 
ducts to convey air flow from the wing or fuselage 
intake to the engine compressor. Compressible, subsonic 
flow conditions usually exist for this application. Ex- 
amples of commercial aircraft with S-ducts include the 
Boeing 727 and Lockheed Tristar. Amongst military air- 
craft, both the General Dynamics F- 16 and McDonnell- 
Douglas F-18 use S-shaped ducts. 

Well designed diffusing S-ducts should efficiently 
decelerate the incoming flow. To achieve appropriate en- 
gine performance, the S-duct must also incur minimal to- 
tal pressure losses and deliver nearly uniform flow with 
small transverse velocity components at the engine com- 
pressor entrance. Unfortunately, the centerline curvature 
associated with S-ducts generates significant secondary 
flows. In addition, a cross-sectional area increase makes 
flow separation likely. Secondary flows and flow separa- 
tion can result in unacceptable duct performance. 

A few investigators have experimentally studied S- 
duct flows, obtaining measurements with pressure probes 
and laser anemometry. 1 ^ Considerable computational 
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effort has also been focused on the complicated three- 
dimensional flows through S-ducts. Researchers at the 
NASA Lewis Research Center recently turned their at- 
tention to complete Navier-Stokes analyses, 5-7 following 
much effort using reduced Navier-Stokes equations. 8-11 
These studies suggest the formation of a pair of counter- 
rotating vortices located within an S-duct. Streamwise 
flow separation was also a dominant feature for S-ducts 
with large centerline displacements and/or cross-sectional 
area increases. From a design point, the ability to pre- 
dict complicated three dimensional flow phenomena with- 
out extensive experimental testing is desirable. However, 
only a limited amount of benchmark data exists for the 
validation of computational fluid dynamics codes applied 
to compressible S-duct flows. 

The objective of this study was to provide a com- 
prehensive benchmark data set for the compressible flow 
through a representative diffusing S-ducL Also, details 
of the flow separation region and the mechanisms which 
drive this complicated flow phenomenon were investi- 
gated. The detailed aerodynamic data are suitable for 
validating computational fluid dynamics methods. The 
experimental data were acquired along with a concurrent 
computational effort, with considerable beneficial syner- 
gism resulting from the combination. This concurrent 
approach also proved to be successful in a recent study 
of circular-to-rectangular transition duct flows. 12, 13 

Experimental Facility 

Detailed data for the flow through a circular cross 
section diffusing S-duct were acquired using the NASA 
Lewis Research Center Internal Fluid Mechanics Facil- 
ity. Complete details of the facility are described by 
Porno et al} A This facility was designed to support tests 
of a variety of internal flow configurations. For this ex- 
periment, atmospheric test cell air was drawn through a 
large settling chamber, passed through the test section 
and dumped into a large exhaust plenum. A schematic 
of the facility is shown in Fig. 1. 


The settling chamber conditioned the incoming flow 
with a series of wires, meshes and screens. First, a per- 
forated spreader cone mixed the inlet flow. Second, a 
coarse mesh conditioning screen reduced mean flow non- 
uniformities. Finally, a honeycomb-screen combination 
removed large scale turbulence fluctuations. A contrac- 
tion section uniformly accelerated the flow from the set- 
tling chamber. An area contraction ratio of 59 ensured 
a low turbulence intensity flow field at the test region 
entrance. 

The test section consisted of the diffusing S-duct and 
two constant area circular extensions. The first extension 
served as the interface between the contraction exit and 
the S-duct entrance. The second conveyed the flow from 
the S-duct to the exhaust region. Each extension was 
76.2 cm long and of appropriate diameter. 

The exhaust region contained a circular pipe, a mass 
flow plug and a sub-atmospheric plenum. The pipe, 
244 cm long, housed the adjustable mass flow plug and 
ensured no downstream influence of the exhaust plenum 
on the test section. The mass flow plug controlled the 
airflow through the entire facility. Flow was simply 
dumped into the exhaust plenum after passing the mass 
flow plug. 

Diffusing S-Duct Geometry 


The diffusing S-duct was intentionally designed to 
incorporate a number of complex three dimensional flow 
features, including separation, which are associated with 
similar configurations. The geometry of the diffusing S- 
duct is shown in Fig. 2. The duct centerline was defined 
by two planar circular arcs with identical radii, R y of 
102.1 cm and subtended angles, 6 max / 2, of 30°. The 
centerline coordinates, given by Eq. (1), are indicated by 
the dashed line in Fig. 2. All cross sections perpendicular 
to the centerline were circular. The duct inlet radius, 
ri, was 10.21 cm. The duct exit radius, r 2 , was 12.57 
cm, which produced an area ratio, A 2 /Ai, of 1.52. The 
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variation of the duct radius as a function of the angle 
9 is given by Eq. (2). The duct is larger than, but 
geometrically similar to, the duct tested by Vakili et 
al. A When discussing results, axial position will refer 
to the distance to the cross-stream planes measured along 
the duct centerline, beginning at the start of curvature. 
Positions within any cross-stream plane are specified by 
the polar angle <£, measured from positive z-axis, and the 
radial distance r from the centerline. 


For 0 <9 < 9 max / 2 
x t \ = R sin 9 
yd = 0 

z c i = R cos 9 — R 


For 9 max / 2 < 9 < 9m ax 

x c \ = 2# sin ( 9 max /2 ) - Rsin ( 9 max - 9) 

Vd - 0 

Zd - 2/? COS (#mar/2) — R — R COS ( 9 max - 9) 

( 1 ) 



The duct was milled from two separate blocks of alu- 
minium and had a final tolerance of ±0.0127 cm. After 
milling, the two halves were mated together, machining 
imperfections were removed, and the surface was pol- 
ished. The interior split line, located on the vertical cen- 
terline plane, was flush to touch and invisible to sight. 



Plane E 

Fig. 2 Half shell of the circular diffusing S-duct. 


Instrumentation and Measurement Techniques 


The flow field was investigated by a variety of meth- 
ods. The internal flow field was measured with calibrated 
three-and five-hole probes. Data were accumulated in five 
measurement planes, each perpendicular to the duct cen- 
terline. The measurement plane locations are presented 
in Table 1. These planes are also depicted by the shaded 
cross sections in Fig. 2. At each plane the S-duct con- 
tained port holes, allowing probe access to the flow. The 
number of port holes and measurement positions con- 
tained in each plane are also listed in Table 1. Reported 
internal flow measurements were concentrated in only one 
symmetric half of the duct. A three-hole probe was used 
at Plane A, since a one-directional velocity field existed 
there. A five-hole probe was used at the other four planes, 
enabling total pressure, static pressure and three compo- 
nents of velocity to be measured. Both yaw-nulling and 
non-yaw-nulling methods were employed. Flow phenom- 
ena were also studied with surface static pressure data. A 
total of 220 static pressure taps were located on the sur- 
face of the duct. The taps were positioned axially at 
constant angles </> = 10°, 90°, and 170° and also circum- 
ferentially at Planes A, B, C and D. All reported aero- 
dynamic data was taken as an arithmetic average of a 
number of measurements. The characteristics of the near 
surface flow were visualized by applying fluorescent oil 
to the surface. For this technique, a powder fluorescent 
dye was mixed with 140 wt. oil and then applied to the 
surface in small dots. Surface streakline patterns devel- 
oped while air was drawn through the duct The duration 
of the test, at the desired flow rate, was ten minutes. 


Table 1 Measurement plane information. 


Plane 

A 

B 

C 

D 

E 

s/d 

-0.50 

0.96 

2.97 

4.01 

5.73 

Ports 

10 

11 

11 

11 

19 

Points 

590 

462 

506 

539 

930 


All pressure data were obtained with a PSI measure- 
ment system. The fully integrated instrument consisted of 
electronically scanned pressure transducers and a micro- 
computer based data acquisition system. Individual pres- 
sure transducers provided high data acquisition rates for 
multiple pressure measurements. The accuracy of each 
transducer was maintained by frequent on line three-point 
calibrations. After acquisition, information was carried 
by Escort, a data routing network, to a storage area for 
later post-processing on Sun SPARCstations. 

An L. C. Smith actuator was used to position all 
probes. Both translational and rotational motions were 
possible. Positions were measured by voltage poten- 
tiometers to ensure accurate location documentation. 
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Inlet Flow Conditions 


Vo Pci 


Inlet flow conditions were obtained from a survey 
of the flow field one radius upstream of the first bend 
(Plane A). All reported tests were conducted with an inlet 
centerline Mach number of 0.6. The Reynolds number, 
based on the inlet diameter and centerline velocity, was 
2.6 x 10 6 . A thin turbulent inlet boundary layer existed. 
A representative inlet velocity profile, plotted in nondi- 
mensional law-of-the-wall coordinates, is shown in Fig. 
3. Comparisons indicated little deviation from a conven- 
tional turbulent profile. Boundary layer parameters are 
listed in Table 2. The boundary layer thickness, 6 , corre- 
sponds to 95% of the free stream velocity. The displace- 
ment thickness, momentum thickness, 6 2 , and shape 
factor, H, were calculated by numerically integrating the 
survey data, implementing a trapezoidal approximation. 


Table 2 Inlet boundary layer parameters. 


Profile Parameter 

Calculated Value 

(6/r i) x 100 

6.95 

(£ i/rj) x 100 

1.46 

( 62 A 1 ) x 100 

1.06 

H 

1.38 


Results and Discussion 

All data presented in this paper are in nondimen- 
sional form. Aerodynamic results represent mean values. 
Pressures are presented as total and static pressure coeffi- 
cients, which are defined by Eqs. (3) and (4). The pres- 
sures p 0 and p represent local values of total and static 
pressure. Inlet centerline conditions define the reference 
states po t ci and p c /. Three dimensional velocity compo- 
nents were converted to local Mach vectors and normal- 
ized by the inlet centerline Mach number, as shown in 
Eq. (5). 



Fig. 3 Typical law-of-the-wall velocity profile. 


C P0 — 


Po,cl Pci 


(3) 


, _ P-Pcl 
Po,cl ~ Pci 


(4) 


M = 


V/c 

Md 


(5) 


For this section the nondimensional pressure coeffi- 
cients will simply be referred to as total and static pres- 
sures. The velocity vector will be presented as normal 
and transverse components. The normal component is 
perpendicular to the measurement plane while the trans- 
verse component is parallel to the measurement plane. 
Since data were taken in only one half of the duct, all 
internal plane results have been mirror imaged for pre- 
sentation. 


Flow Visualization 

Major flow field characteristics were visualized by 
streaklines of fluorescent oil on the duct surface. Three 
conclusions were made from the streakline patterns. First, 
the flow was indeed symmetric. This was ascertained 
by applying different color dye to each symmetric duct 
half. The dye from each half remained separate. Second, 
a large region of separated flow existed, as strikingly 
indicated in Figs. 4 and 5. The entire duct section is 
shown in Fig. 4, while Fig. 5 is a close up of the 
separated region. The free stream flow is from left to 
right and each photo shows only one symmetric half of 
the observed streakline pattern. The separated region 
consisted of two saddle points, occurring on the duct 
split line (<£ = 180°), and two spiral nodes, lying in 
each symmetric duct half. The entire separated region 
was located on the lower wall. The onset of separation 
(upstream saddle point) was located at s/d = 2 . 02 , while 
reattachment (downstream saddle point) occurred at s/d = 
4.13. Third, boundary layer cross flows were present. In 
the first bend, streaklines were driven towards the lower 
surface (</> = 180°). In the second bend, the streaklines 
near <j> = 90° diverged. The upper streaklines converged 
towards the top of the duct, while the lower streaklines 
continued towards the lower split line (</> = 180°). 

A qualitative assessment of the flow in the plane of 
symmetry (x-z plane) was acquired by temporarily intro- 
ducing a splitter plate, which divided the two symmetric 
duct halves. The plate was installed for only one test 
in order to produce streakline patterns. The S-duct sur- 
face results given above were obtained without the splitter 
plate in place. The splitter plane results are assumed to 
be qualitatively correct, even though the presence of the 
flat plate introduced additional shear into the flow. The 
streakline patterns produced on the splitter plate are pre- 
sented in Fig. 6 . Free stream flow is from left to right. 
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Fig. 4 Surface oil flow patterns on one symmetric half of the duct surface. 



Fig. 5 Close up of the surface oil flow' patterns developed in the separation region. 
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Fig. 6 Surface oil flow patterns on the centerline splitter plate. 


The location of separation and reattachment are identifi- 
able. The streaklines follow the duct curvature in the first 
bend. However, blockage caused by separation forced 
the streaklines to deviate from the curvature of the second 
bend. The reversed flow region began thin and appeared 
to grow until reattachment. A strong vortical structure 
followed reattachment. Viscous effects associated with 
the splitter plate can be seen by the extensive bending 
of the streaklines at the exit of the duct. Cross-stream 
pressure gradients helped drive the boundary layer fluid 
on the splitter plate towards the top of the duct. This 



was most noticeable at the exit where the splitter plate 
boundary layer was the thickest. 

Surface Static Pressures 

Symbols in Fig. 7 represent the static pressure vari- 
ation with axial distance for three circumferential loca- 
tions. The region of streamwise separated flow, deduced 
from flow visualization, is also shown. The effects of 
streamline curvature and diffusion are clearly indicated 
during the first bend by the pressure difference between 



Fig. 7 Axial distribution of surface static 
pressures for three circumferential positions. 
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Fig. 8 Circumferential distribution of surface 
static pressures at four axial locations. 





Plane B 
Plane C 
Plane D 


s / d = 0.96 
s/d = 2.97 
s / d = 4.01 



Plane E 


Fig. 9 Static pressures in three interior planes. 


<j) = 10° and 170° and overall pressure rise, respectively. 
The influence of separation is shown by the constant val- 
ues of the (f> = 90° and 170° static pressures between 
2.0 < s/d < 3.2. The proximity of these two curves, 
between 2.0 < s/d < 3.2, indicate minimal change in 
flow speed or direction in the lower half of the duct. The 
blockage, caused by the separated region, also helped to 
produce a favorable pressure gradient at <j> = 10°. Static 
pressures at the three angles are nearly equal at s/d = 2.6, 
suggesting a constant static pressure distribution through- 
out the cross-stream plane. The static pressure rose again 
for s/d > 3.2, even though the flow remained separated 
beyond this region, because the duct continued to diffuse. 
After reattachment the pressures distributions converge, 
although not plotted, to Cp = 0.466 far downstream 
(s/d « 9.0). 

The circumferential distribution of surface static 
pressure in Planes A through D is presented in Fig. 8. 
The pressures for Plane A are nearly equal indicating no 
influence of the duct curvature at this upstream section. 
The data for Planes B, C and D all reflect the presence of 
streamline curvature. The pressures measured in Plane B 
are maximum at <\> = 10° and continually descend until 
(j) = 170°. The values are negative for 0 > 110°. Sepa- 


ration has not affected the flow field yet. In Planes C and 
D the cross-stream pressure gradient has reversed. How- 
ever, the influence of separation can be seen by the slight 
drop in pressure at large values of <j>. The surface static 
pressures also indicate the separation region has reduced 
in size from Plane C to D. 


Internal Aerodynamic Data 


Contours of static pressure are shown in Fig. 9 for 
the Planes B, C and D. Static pressure data indicated a 
nearly uniform cross-stream static pressure distribution 
at Planes A and E and therefore are not presented. The 
values of static pressure very near the duct surface in Fig. 
9 are in excellent agreement with the data obtained from 
surface static pressure taps shown in Fig. 8. 

Total pressure distributions for the five planes are 
presented as contours in Fig. 10. Expanded views of the 
distributions of total pressure for Planes A and B are also 
included in order to help distinguish the boundary layer 
contours. The convection of boundary layer fluid away 
from the surface into the core flow is quite pronounced 
at Planes C, D and E. 
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Plane E 


Fig. 10 Total pressures in the five measurement planes. 


Contours of the normal components of velocity at 
all five planes are displayed in Fig. 11. An expanded 
view of the inlet boundary layer is presented for Plane 
A. A large region of low momentum fluid can be seen 
developing in the last three planes. Planes C and D lie 
within the region of separation. However, at the lowest 
contour level of the normal velocity component it was 
impossible to ascertain with a five-hole probe whether 
the flow was reversed. The distributions of the normal 
velocity component and total pressure are quite similar. 

Finally, the transverse velocity components are plot- 
ted in Fig. 12 for the last four planes only. Data for the 
inlet plane (A) is not shown since cross flows were not 
detected there. The generation of strong transverse flows 
within the duct can be easily seen. 

The flow field at Plane A corresponds to a develop- 
ing pipe flow. The results in Fig. 3 verify the boundary 
layer is fully turbulent. Both the normal velocity and to- 
tal pressure contours, shown for Plane A in Figs. 10 and 
11, show no circumferential variance. 

As the flow progressed downstream, static pressure 
distributions developed from core flow adjusting to duct 
geometry and interacting with the separated region. Ini- 
tially the S-duct deflected the incoming flow downward, 


which created large cross-stream pressure gradients. The 
static pressure distribution at Plane B in Fig. 9 shows 
the maximum static pressure value at the top of the duct 
(<!> — 0°) and the minimum at the bottom of the duct 
(</> = 180°). The normal velocity distribution at Plane 
B in Fig. 1 1 again shows the effect of streamline cur- 
vature. At the top of the duct, where the cross-stream 
pressure gradient is the greatest, there existed an adverse 
streamwise pressure gradient that decelerated the flow. 
Likewise at the bottom of the duct, where the minimum 
cross-stream pressure is located, there existed a favorable 
streamwise pressure gradient upstream which accelerated 
the flow. This can be seen in Fig. 11, where the lo- 
cal Mach number exceeded the inlet Mach number near 
the bottom of the duct. The transverse velocity compo- 
nents in Fig. 12 show that the cross flow occurred only 
near the duct surface. This was caused by the pressure 
gradient turning the lower momentum fluid towards the 
bottom of the ducL This was also detected in flow vi- 
sualization. There was little change in the total pressure 
distribution at Plane B, given in Fig. 10, other than a 
slight thickening of the boundary layer. The boundary 
layer thickness at Plane B did also vary slightly with cir- 
cumferential position. The boundary layer thickness was 
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Plane A 

s / d = 

-0.50 

Plane B 

s / d = 

0.96 

Plane C 

s / d = 

2.97 

Plane D 

s/ d = 

4.01 

Plane E 

s/ d = 

5.73 



Plane E 


Fig. 11 Normal component of Mach vector in the five measurement planes. 


greatest at </> = 0° and least at <j> — 180°. These phenom- 
ena can be attributed to the streamwise pressure gradients 
which accelerated the flow near the bottom of the duct 
and decelerated the flow near the top of the duct. 

Near the middle of the S-duct the static pressure dis- 
tribution should become nominally flat as the centerline 
curvature changes. In the second bend the orientation 
of the cross-stream static pressure distribution should re- 
verse, so that the lowest static pressure is near the top and 
the highest static pressure is near the bottom of the duct. 
Both Planes C and D in Fig. 9 show this reversal. The 
static pressure distribution at Plane C, which was located 
closer to the middle of the duct and in the separated flow 
region, was flatter than the distribution at Plane D. Both 
distributions indicate that an area of minimal change in 
flow speed or direction was present in the lower half of 
the duct, confirming the previous conclusion speculated 
from the surface pressure data. Normal velocity compo- 
nents (Fig. 11) and total pressure distributions (Fig. 10) 
show that a low velocity region existed in Planes C and 
D in the lower half of the duct. This region was cause by 
reversed flow and the convection of boundary layer fluid 
toward the lower surface. Strong pressure-driven cross 
flows existed near the lower duct surface in Planes C and 


D (Fig. 12). These cross flows continually drove bound- 
ary layer fluid into the separated region. Enough low 
momentum fluid was convected upwards to form a lo- 
cal minimum of total pressure and normal velocity away 
from the surface. This suggests that a vortical structure 
was developing between Planes C and D. A significant 
amount of flow deflection at the centerline in Plane C 
is shown in Fig. 12. The deflection was a direct con- 
sequence of the blockage created by separation. This 
confirms the results obtained from the flow visualization 
on the centerline splitter plane. The core flow returned 
to the nominal streamwise direction by Plane D, which is 
contrary to the splitter plane patterns. This inconsistency 
is a result of the pressure gradients which drove the build 
up of boundary layer fluid on the splitter plate toward the 
top of the duct, as previously explained. 

By Plane E the free stream flow returned nominally 
to the x-direction. Cross-stream static pressure gradients 
were eliminated. However, before reaching the exit plane 
(E), static pressure gradients were strong enough to re- 
verse the direction of the boundary layer fluid on the 
upper surface of the duct, as seen in Plane E (Fig. 12). 
A large pair of counter-rotating vortices, located in the 
lower half of Plane E, also evolved. These vortices con- 
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Plane A 


Plane B 

s/ d = 

0.96 

Plane C 

s/d = 

2.97 

Plane D 

sl d = 

4.01 

Plane E 

s/ d = 

5.73 



Plane D 


Plane E 


Plane E 


Fig. 12 Transverse components of Mach vector in four measurement planes. 


tinually convected the low momentum fluid of the bound- 
ary layer towards the center of the duct. This convection 
degraded both the uniformity and magnitude of the total 
pressure profile. The region of low momentum fluid in 
Plane E (Fig. 10 and 11) extend above the center line. 

Summary 

The flow through the S-duct evolved from a strong 
interaction between the boundary layer and the core flow. 
Results show the flow was symmetric about the x-z plane. 
A large region of streamwise flow separation occurred 
within the duct Duct curvature induced strong pressure- 
driven secondary flows, which evolved into a large pair 
of counter-rotating vortices. These vortices convected 
the low momentum fluid of the boundary layer towards 
the center of the duct, degrading both the uniformity and 
magnitude of the total pressure profile. The benchmark 
data collected should be beneficial to inlet designers and 
is suitable for code validation. 
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